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ABSTRACT 


In a welded structure, shrinkage and externaleconsecoueie 
are closely related. Shrinkage decreases as Degree of Con- 
straint increases. 


In the first part, this paper presents: 


(1) a formula that can be used to find the transverse shrin- 
kage of a free joint in the case of a butt weld, 
(2) the definition of the degree Of Gonstrainte Gferes = eames 
i structures, 
(3) the empirical correlation between transverse shrinkage 


and degree of constraint. 


In the second part, this paper presents: 

(1) analytical formulae for the deduee otmconsteatmminot ms ene 
Simple Joint contvcuratloncr 

(2) the numerical method for the degree of constraint of 
other joint configurations and thewtabulatvoenm.s. sen omeoe 
Sults obtained. 


In the last part, the use of the degree ©f const uaame 
to determine the cracking susceptibility of a weldecmemiane= 
ture 1s suggested and a method of experimentation to verify 
the numerical results obtained in the second part is proposed. 
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SUMMARY 


Distortions and residual stresses unavoidably accompany 
every welded structure. They are not independent of each 
other. In the case of an Unrestrained Structure, ee ono an. 
kage distortion is determined by the welding heat and can be 
determined. However, practical structures are seldom com- 
pletely free and joleee 1s performed more or less under ex- 
ternal constraint, therefore it seems that if factors rela- 
ting to welding heat are definite, the ratio of shrinkage 
distortion under external constraint to the one ina free 
welded joint 1s uniquely decided by the degree of external 


constraint to which weldments are subjected, or: 


(shrinkage under external constraint, Si) = (shrinkage 


in free welded joint, Sy) xX ee 


where ee is a function depending on the degree of constraint. 


Sie is the transverse shrinkage of the unrestricted joint or 


free JOInt. Fouvea Duties euieeees has the following expres- 


GE 


sion: 


where: 


by 
I 


sectional area of the groove Of Ene buee jomme- 


=i 
Il 


plate thickness, 


= 
II 


weight of deposited metal per unit weld length, 








W, = weight of deposited metal per unit weld length per 
welding of each pass, and 


C, and C. are coefficients depending on arc voltage, arc 


aC 2 


current, arc efficiency, filler metal diameter, and 


melted weight of rod per unit current and unit time. 


For example with the IJlumenite type filler rod of 3.2 mm 
diameter, with weld current of 120A and welding speed of 0.3 cm/s, 
it has been found that in the C.5.G. unit system: 


Cy = 0.0960, C. = 0.0416 


In this order of idea, Masubuchi has defined the degree of 


constraint as: 


Uniform transverse stress, oy 
@) 
Degree of Constraint,K = ———————__ (1) 


Average Transverse Displacement [v], 


where Iv], is the average of the transverse displacement taken a- 
long the weld length &, and has found a graphical correlation 
between transverse shrinkage and degree of constraint K (refer- 


ences 1,2). 


Later, Watanabe and Satoh, after analyzing the experimental 
results obtained by Masubuchi and others, have found a correlation 


of the form: 


ep) 
Gi 


7 1 
= ee (25) 


te 1+ 0.086K2° 


therefore, in order to determine the shrinkage distortion, the 


Cp) 


value of K must be known. 


On the other hand external Gonstraint promotes sires aa 


fact, the level of stress can be found by combining (| Syaeae (2): 








oo = KI 
Yo y 
and for a butt weld, S, = iv], for an actual weld. Therefore: 
K O is the average reaction 
O - Sig me Yo stress along the weld due 
Yo 1+0.086K~° to constraint 


If the external constraint is too high, residual stress can reach 
a dangerous level where cracking can occur. Here again, it is es- 


Sent lal Acour pow. 


In a few cases, K can be found analytically. But in most 


practical cases, it has to be computed numerically. 


The object of this study is to compute numerically the values 
of K for some test specimens and to propose a method of experimen- 
tation to verify the computed results. Also the degree of con- 


straint in the case of patch weld has been determined analytically. 


The numerical method used is the finite element method in 
which the plate structure is divided into small elements in each 
of which static equilibrium and geometrical compatibility have to 
be satisfied. In this case, for each specimen the values of K 
are computed as a function of R = - where £ is the weld length 
(in this case, it is the length upon which the uniform stress o ; 
acts) and L is the length of the slit. The results are presented 
EID (E is 


the modulus of elasticity) curves versus R = £/L in the attached 


as the non-dimensionalized degree of constraint K = 








diagram (Figure 9). 


One exception is made in the case of the Lenigh teste spe- 


cimen where the weld is made all along the length of the slit 


| 
| 





SS a RI - —— 


without saw cut with saw cut 


(2=L, R=1.0). In these cases, experimental data are available 


for comparison: 


a eee a A A i A 


: | Kcomputed “measured 


K * K 2. 
computed measured K computed 
No Saweut 45.3 Kg/mm*-mm 44 Kg/mm*-mm 22 8m, 
| 
With Sawcut (28.9 Kg/mm*-mm /27 Kg/min? rn 6.70 


ccc rence neces arent ethnics engine cpp emer ammenities 
* 
Assuming E = 19.9 x 10? Kg/mm?* 


The object of the Lehigh test specimen is to find the critical 
value of K or range of values of the degree of constraint above 
which weld cracking may develop, below which it may not. The ana- 
lytical results obtained in the case of a patch weld ,which consists 
of a circular disc welded to a plate, are fon the easeuoreaee 


with large dimensions (infinite plate): 





Patch Weld 


If the weld can be made along the whole circle in one pass, 


then: 


R is the radius of the plate, 
Oo. is the radial stress of the perimeter of the disc, and 


1s the radial displacement at the weld. 


5 (2) If the weld can be made within a sector 2a, then: 





— Zou 2-E (Ae) 


RJ 20 


5) where (hall 5, is the average radial displacement over the welded sec- 


tor 20 and £(20) 1s a series function of 20 with £(0) =n) one 


. dk 
f) the series converges as —. 
n 


Because of the linearity of the relations in elastic defor- 
mation, the displacement 1s inversely proportional to the modulus 
of elasticity E, and K, which is defined as a stress divided by 
the corresponding average displacement, iS proportional to E. 
Therefore a new factor, called specific degree of constraint K', 


and defined as: 


Ki = 


A 


is independent of the linear elastic material and dependent only 


on the geometry of the weld. 


Therefore experiments to determine the values of K to verify 


the computed results can be conducted on any linear elastic ma- 





terials, preferably the ones that have small E and are not expen- 
Sive. Plastic materials are proposed, and a mechanical way to 


stress the specimens is to drive two wedges to make an assembly 


c7¢ ae wedge shape 


with parallel faces into the slit. Each specimen will have 


several sets of wedges with different widths & so that K' can be 


£, 


f) found as function of R = rc Strain gauges put near the edges of 








the slit and visible marks made at the immediate vicinity of the 
edges for optical measurements by microscope or comparator are 
used to determine stresses and displacements. In these experi- 
ments the material at the loading zone is in compression instead 
of tension as in the welding. It is of no consequence as far as 
the material is linear elastic in compression. it 1S eacruaielne 


and the value of E is between 3 to 5 x 10° psi (in compression). 


The problem of shrinkage control could be conceived as a 
problem of choice between two alternatives, based on the degree 


Ot (COMmS tar aime ee 








O 
Yo 
K = on 
[vl], 
and the Watanabe-Satoh relation: 
ne 
Fer Ste 140.086 K%8? 


Sif is the transverse shrinkage of a free joint and for a butt 


weld, §S Iv], at the weld. These alternatives are either to 


‘cen 
limit the shrinkage by increasing the external constraint, and 
thus accept a higher value of K and hence a higher level of re- 
Sidual stress or to reduce the residual stress level by relaxing 


the external constraint thus reducing the value of K and to accept 


a larger value of shrinkage distomelont 


In the problem of cracking control, it is essential that the 
external constraint will result in a value of K below the range 
of critical degree of constraint Kas characteristic of each ma- 


terial and each weld type. 





In these two problems, the important factor is the degree 


of constraint K. 








INTRODUCTION 
I. General 


Welding is an efficient method of fabricating ss eimte.. 
ture. Various products ranging from miniature electronic 
components to bridges, ship hulls, rocket motor cases have 


been made by welding. 


One of the troublesome problems that accompany the con- 
struction of welded structures is shrinkage distortion. The 
more complex the structure is, the more involved the problem 
becomes. Shrinkage distortion can cause mismatch of joints 


which leads to the possibility of welding defects. 


The correction of weld distortion is costly and in 
some cases impossible. It is therefore desirable to develop 
some techniques to predict somehow the approximate values 
of the shrinkage distortion in order that palliative mea- © 
sures can be devised that will neutralize or reduce the 


effect of shrinkage distortion. 


Another problem associated with welding is cracking. 
Cracks may form as a result of the welding operation and can 
occur within the weld metal or the base metal in the heat- 


affected zone. 


There are two types of cracking that are experienced in 
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welding: hot cracking and cold cracking. Hot cracking is be- 
lieved to take place during the solidification of the weld- 
ment and can occur in the weld metal or in the heat affected 
zone. Hot cracks are of intergranular type. Cold cracking 
occurs at much lower temperatures than does hot cracking, 

in steels at temperatures below the start of the austenite- 
martensite transforiuation. Cold cracks may be formed during 
or after the welding operation; cold cracks are of trans- 
granular type. Cracking, hot or cold, depends on many fac- 
tors such as chemical composition of the base and filler 
metals, mechanical restraints, welding condition (heat,input). 
For given welding conditions, chemical compositions of the 
base and filler metals, it is desirable to develop a way of 
determining quantitatively the degree of restraint that may 


PrOMOue se sacle 


it. Technical Backonoumel 


Distortion in a welded structure may be determined as 
eeeEunc Clow et structerel parameters, material ee ernctare and 
fabrication paramerers, 

The structural parameters include the geometry of the 
structure, the shape of the joining boundary and the type of 
the joint. 

The material parameters are the nature of the base and 
filler materials. 


The fabrication parameters include the welding process, 


the heating procedure, the welding sequence, the degree of 









IE L 


constraint. 

The degree of constraint, which will be the principal 
object of this study, is a parameter which characterizes the 
external mechanical constraint to which a weld is subjected. 
For example, the degree of constraint of a weld made through- 


out along the median line of a rectangular plate, case (a), 
——— 





ee es a et ta cm 


is different from that of a weld made only partially along 
this line, case (b). In fact, in case (a), the weld can be 
considered to be free while in case (b), it is not. And in 
case (a) the shrinkage is larger than in case (b) and we can 
imagine that there should be a relation between shrinkage and 
degree of constraint. 

According to Masubuchi (reference 2), to analyze weld 
distortion, it 1S necessary to establish analytical relation- 
ships among these three sets of parameters and distortion. 

For a Simple butt weld as is the case of this study, the 
dimensional changes produced in the structure by each weld are 
to be determined. This can be done by: 

(1) Analyzing the heat flow, 


(2) Analyzing the thermal stresses during welding to 


determine incompatible strains, and 


(3) Determining the dimensional changes. 
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In ‘Fusion welding, a weldment is locally heated by the 
welding heat sources. During the thermal cycle, the weldment 
is subjected to thermal stresses. When the weld is completed, 
incompatible strains are created in the region near the weld. 
Incompatible strains, including dimensional changes associated 
with solidification of the weld metal, metallurgical transfor- 
mations, and plastic deformations, are the sources of residual 
stresses and distortion. When welding processes and parameters 
are changed, the heat flow pattern is also changed causing a 
change in the distortion of incompatible strains, hence in 
shrinkage and distortion. 

The problem of determining the distribution of incompatible 
strains is extremely difficult. When a material undergoes plas- 
tic deformation, the stress strain relationship is not linear. 
Furthermore, plastic properties of the material change with 
temperacure. 

When the incompatible strains are known,theoretically or ex- 
perimentally, the problem of determination of dimensional changes 
can be handled analytically. Moriguchi has developed a funda- 
mental theory of stresses caused by incompatible strains, and 
Masubuchi has applied Moriguchi's theory to the study of resi- 
dual stresses and distortion due to welding. 

Assuming that the dimensional changes in welds are determined, 
either analytically or experimentally, the next step is to deter- 
mine the distortion induced in the structure by these dimensional 


changes. Although plastic deformation is produced in small areas 
| 














ples 


near the weld, most of the remaining material in the structure 
is elastic. Consequently, the induced distortion can be ana- 
lyzed by using the elastic assumption. Hence the relations used 
to determine the induced distortion are independent of fabrica- 
tion parameters and depend only on well-established material 


parameters. 


A. Changes of Temperature and Stress During Welding (from refer- 
ence 2) 


Figure 1 shows schematically how residual stresses are formed 
in a weld. Figure la shows a bead-on-plate weld in which a weld 
bead is being laid at a speed v. O-xy is the coordinate axis; 
the origin, O, 1s on the surface underneath the welding arc, and 
the x direction lies in the direcrion o2ave dina 

Figure 1 shows temperature distribution along several cross 
sections. Along Section A-A, which is ahead of the welding arc, 
the temperature change due to welding, AT, is almost zero (Figure 
lb-1). Along Section B-B, which crosses the welding arc, the 
temperature et is very steep (Figure WD) Along 
Section C-C, which is some distance behind the welding arc, the 
distribution of temperature change is as shown in Figure 1b-3. 
Along Section D-D, which is very far from the welding arc, the 
temperature change due to welding again diminishes (Figure l1b-4). 

Figure lc shows the distribution of stresses along these 


Sections in thevxsdireceion, OL: Stress in the yeaqurec ton: Gt 


and shearing stress, Tay! also exist in a two-dimensional stress 
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NiO Stress XQ 
ase: eae 
D * 
"1 Section A-A 
| Ls Cot ate teteiatete A . 
U | ; . 
4 Y : | ° 
Melted region | 
Region ae | : 





where plastic : | 
deformation : 2. Seclion B-B . 
occurs | 

during welding 





3. Seclion G26 
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: : : 7 ij! 
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eS 


* 
field (Figure la). 


Along Section A-A, thermal Stresses due to welding are al- 
most zero (Figure le-1). The stress distribution alongeseeerem 
B-B is shown in Figure lc-2. Stresses in areas underneath the 
welding arc are close to zero, because molten metal does not 
Support loads. Stresses in areas somewhat away from the arc are 
compressive, because the expansion of these areas is restrained 
by surrounding areas that are heated to lower temperatures. Since 
the temperatures of these areas are quite high and the yield 
strength of the material is low, stresses in these areas are as 
high as the yield strength of the material at corresponding tem- 
peratures. The amount of compressive stress increases with in- 
creasing distance from the weld or with decreasing temperature. 
However, stresses in areas away from the weld are tensile and 
balance with compressive stresses in areas near the weld. In 


other words, 


fo, - dy = 0 (a) 


kx | 
aCKOSS SeCCLIGne EL. Thus, the stress distribution along Section 


BB is as shown in Figure lc-2. 


Stresses are distributed along Section C-C as shown in 


* 
In a general three-dimensional stress field, six stress compo- 
Nents ,~ 0 Gy Gen 6 ee ee @xXLst. 
mc y Z. XY ZN, ZK 
k* ; ; } 
Equation (a) neglects the effect of Ss and Ty on the egquili- 


DEM scOne i tile@me 
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Figure lc-3. Since the weld-metal and base-metal regions near 
the weld have cooled, they try to shrink causing tensile stresses 
in areas close to the weld. As the distance from the weld in- 
creases, the stresses first change to compressive and then be- 
come tensile. 

Figure lc-4 shows the stress distribution along Section D-D. 
High tensile stresses are produced in areas near the weld, while 
compresSive stresses are produced in areas away from the weld. 
The distribution of residual stresses that remain after welding 
is completed are shown in the figure. 

The cross-hatched area, MM', in Figure la shows the region 
where plastic deformation occurs during the welding thermal cycle. 
The cross-hatched area near the origin 0 indicates the region 
where the metal is melted. The region outside the cross-hatched 
area remains elastic during the entire welding thermal cycle. 

Because of the difficulty in determining the distribution 
of incompatible strains, no analysis has yet been developed to 
trace the com of two-dimensional thermal aimessos Gumine 
welding and to determine distributions of three residual-stress 
components, Ow Oy" and T yy" In other words, no analysis has 
been made in which both heat flow and Stress ivelas=a-c meme cecs 
aS two-dimensional problems. In all studies conducted so far, 


the problem has been simplified in some way. 


B. Shrinkage Distortion 


Shrinkage distortion in welding can be considered as the 











IEG 


& result of the combination of tworqnoups er sracrers,. 


(1) Factors related to welding ware, ane VoOlEage yeaa 


current, welding speed, type and size of electrode, 


(2) Factors related to external constcaine. 


os 


lj= BErect Of Weldinean lect 


The effect of welding heat on the transverse shrinkage in 


the case of butt weld free of external constraint is known in 





lave AeengNt= 
Sie = C1 52 Log a ie Qa) Ae 

where 

Sif is the transverse shrinkage of the free joint, 

A is the sectional area of groove of the butt joint, 

b 1s the thickness of the plate, 

W is the weight of deposited metal per unit length of 

the weld, 


W,. 1s the weight of deposited metal per unit length Oe 
the weld in one pass, and 


C, and C 


Ik 


2 


are constants depending on the arc voltage, arc 


intensity, heat efficiency of the welding arc, size and 


type of the filler metal and the melted weight of the 


rod per unit current and unit time. 


2 -Effect of External Constraint 


In practical work, welding is performed more or less under 


external constraint. In this case, shrinkage distortion in weld 


metal is depending also on the mechanical constraints and smaller 


/ 








rae than in free joints. In their work, Kihara and Masubuchi at- 
tempted to determine quantitatively the degree of constraint and, 
in a slit weld, to find a relation between transverse shrinkage 


and degree of constraint. 


QO is the center 
of the plate 





Accordingly, the degree of constraint K is defined as: 
K = —t2 | (2) 
[vl], : 


where £ is the length of the weld in a slit of length L, i is 
Oo 


the uniform stress applied along the weld length 2& and alr is 





the corresponding mean value of transverse displacement over the 
portion of the slit where the load is applied. The physical 


meaning of K is that when uniform transverse stress oO is ap- 
Oo 


plied along the part of the slit between X=Xy and X=Xo¢ 


= (]x5-x, | = 2) displacement v will occur along the slit. The mean 


value of the transverse displacement [v], defined by: 


x 
2 
_ il 
Iv], = | Vax is related to ce Dye 
x1 * ¥ . 
= K 
Ie Iv] 


In the case of a slit in an infinite plate, and using the analogy 
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® between the residual stress in this case and the vortex theory 


in Hydrodynamics, they have shown that (see references.1 and 2): 





Wi eal 
NS 0) ae (3) 
with 
07 O> 
F= 5 [ ff sinOsinmdQ] ” 
n=l - Oy 
where 
QL = X5 ~ Xy5 57%] 
me 
©, = cos es 
ty Z 
x 
oe = COS : =>) 
L/2 
The value of K as a function of R= = has been computed and the 
= K TT R 
curve K = => xX = 


E/ vy, F(R) has been drawn in Figure 2 for the case 
where the weld is symmetric, that 1s x, + xX, = 0; ana 0. = TO, 
the encircled dots on the diagram are the value oL 7h — SE found 


experimentally. We note that K has the dimension of a stress 





divided by a length and therefore K is dimensionless. [In the 
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experiments described in reference 1, the following results have 


been ebrained: 


eee 5 an : ’ } 
Slit ee Length K leh Il | 





Q nen 
| L (mm) | 2 (mm) y (Kg/mm?/mm) | (Kg/mm? /mm) A 
| | 
| 513 | 1g 3 0.337. 51.6 steels 0.905 
340 | 164 (0.488 | 56.4 62.06 0.740 
174 | 174 1 000 76.8 1s 0.636 
523 174 0.330 52.0 40.34 | 0.915 
349 175 10.502 55.2 60.46 0.730 
A 174 1.000% 76.8 LAs OG 0.636 
510 163 0.320 38.3 41.37 0.930 
— 176 176 LOC 76.0 119.89 , 0.636 
fia 171 1.000 VED 123.39 390.636 
340 168 0.494 57.0 62.06 |. 0.735) 
172 172 1.000 77.8 122.67 | 0.636 
| 337 1GD 0.510 55.5 62.61 | 0.725 
| 165 165 1.000 81.0 12 LE 0.636 
/ 202 202 le OOlG 66.2 104.46 | 0.636 
224 224 1.000 59.6 94.20 | 0.636 
129 129 1.000 ee. 5 163.57 0.636 
96 96 1.000 139.1 219.79 0.636 
1 390 149 0.382 45.6 54.10 ore 40 
| 290 150 0.518 Soe 72.76 | 0.720 


In their report, Kihara and Masubuchi did not mention whether the 
welds are symmetric or not in the case of x<l1.00. But it has 
been confirmed! by Masubuchi that the welds are made symmetrically. 


The length of the plate is L, = 1200 mm, its width is B = 800 mn, 


0 


its thickness is h = 19 mm. ‘The Slit Gerais eco nee 


accompanying figure. 
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The experimental values of K are higher than those given by 
the theory, instead of the contrary, since the plates are actually 
finite and the shrinkages are expected to be larger conseguently. 
There are several reasons for this. One as) that heu oueeeetaaees 
be perfectly symmetrical, while the values of K and K are loweer 
in the symmetric case. Another is that the methods of measuring 
the shrinkages as used by the authors give a lower value than 
they are actually, since the width of the slit is about 12.8 mm , 
while the shrinkages are measured along parallel lines which are 
45 mm. apart and the value of the shrinkage at that distance of 
the weld (2mm - 42. 8mm 6.1mm) is expected to be smaller than 
they are at the weld. 

Using this degree of constraint thus defined, Masubuchi has 
found a correlation between the mean shrinkage and Kina slit 
type specimen as reproduced in Figure 3. 

Expanding this idea, Watanabe and Satoh (reference 3) later 
determined the value of K for other configurations such as the 
H-type and the circular-ring type. Then they found an empirical 
correlation between the transverse shrinkage Sy and the degree 


of constraint K for various types of weld in themeorne 


S 1+0.086K0° 8? C 


eats 














being the transverse shrinkage of a free joint which has been 


Ste 
defined in the previous part. And they conclude (reference 3) 


that, "It may be said...that these data are approximately repre— 


sented by one curve despite the fact that various conditions are 


different from one to the other. It can be concluded, therefore, 
S 
that thewaunetren ie (= a) +..0r the ratio of shrinkage dis- 
Gi 


tortion under external “restraint -£o. ene ona] eee 


welded joint 1s decided by K and it is independent of the other 


© 
COnd1ErOns, Bieedmaeneew = = f(K) has been reproduced in 
ete 


Figure 4 with some experimental data for comparison. The value 


of S in the case Of bute we lems: 


ieee 
7 2 
hie Wo h? 


A; sectional area ©f groove or bUGEs jeune 
h: plate thickness 
W: weight of deposited metal per unit length 


Wo: weight, of deposited metal per unit length per 
welding of each pass. 


In the C.G.S. unit system and with the [lmenite type, the co- 


efficients C, and C., have the following values: 


1 2 
Welding Current Welding Speed Size of Electrode Cy C. 
(A) (cm/sec) (diameter in cm) 
120 Os: 32 0.0960 0.0416 
150 Os | 4 0.1021 02 G5e4 
PAY Oe 5 0.1530 0. 07 45 
260 0.3 6 0.1249 0 UGo0 
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FIGURE 3. RELATIONSHIP BETWEEN DEGREE OF CONSTRAINT AND 
| TRANSVERSE SHRINKAGE IN A SLIT-TYPE SPECIMEN 
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The correlation obtained by Watanabe and Satoh presented in 
the previous text enables the prediction of the transverse shrin- 
kage and the level of transverse residual stress of a plate due 
to given welding conditions if the degree of constraint of the 
weld is known, provided that this empirical correlation is cor- 
rect. The Watanabe-Satoh correlation has been obtained from a 
set of experimental data conducted on three types of joints: the 
Straight slit type, the H-slit type and the circular ring type 


Specimen. Some values of K given by Watanabe and Satoh are: 


Se A AS A te A A IG A 


oe | 
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=e 
apes 


| | 
. e— ie —— K == 
“~ | 


(for an infinite plate) 





si "I 
| 
ee on ee E 1 
+ a 
| a 
7 2 
| 
| _— . (No information was given as 
a ann | to the width of the H bran- 
| ! Ls | | ches,) the lengthwot sence 
a a ae | a plate) 
ea 6 - - 


H-Siit Type 
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The straight slit formula has been found by Masubuchi in refer- 
ence (1) or (2). The circular ring formula can be found by using 
the elastic theory (Reference: Theory of Pelastrecrt ey oye eee sie micer 
2nd edition, problem No. 4, page 126). The H-slit formula can be 
found by uSing the elastic theory, according to Watanabe and 


Satoh. 


In these three cases, the geometrical boundaries are rela- 
tively simple and analytical formulae are possible. Unfortunately, 
with the probable exception of the ring type, these specimens are 
not easy to make. Usually to make a slit, it would be easier to 
begin by drilling two relatively large circular holes at the ends 


and then cutting the plate along a line joining their centers. 


Also, the dimensions of the plate are finite. Then the 
region is not single-connected any more and analytirca seo ume 
tion becomes very difficult if not impossible, and numerical 


methods should be used. The same remark can be applied to actual 


Structunesc: 
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The purpose of this study is to compute the value of the 
degree of constraint K for some types of specimens, to propose 
a method of experimentation which, when performed, will provide 
a means of verifying the numerical method used in the computa- 


ELON CO Lek 


This numerical method of computing K provides also a method 
of quantitatively evaluating the cracking susceptibility of a 
weld. In this case, some experimentation methods have been pro- 
posed (see reference 4, page 32-39), such as the Lehigh test, in 
which the values of K are measured for different specimens. If 
the value of K is too high the plate will crack and there should 
be a critical value K., Or a range of values, of the degree of 


constraint K for which: 


K below the critical range, the plate will not crack, and 


K above the critical range, the plate will crack. 


The critical range of K is determined experimentally and 
could be consicered as a material and weld type property. Then 
if a weld is to be made for a given plate material, its ability 
to crack can be predicted by computing the value of K of the 
weld to be made and thus its feasibility can be decided quanti- 


tatively before the weld is actually made. 
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PROCEDURE 


I Mathematical Approach 


In most cases, a numerical method will be used which is 
based on the finite element method of structural analysis de- 
veloped recently and a summary of which can be found in reference 
(5). The computer programs for the finite element methods in 
structure analysis are made available by the Civil Engineering 
Department of M.I.T., whose user's instructions can be found in 


reference (7). 


The basis of the method is to divide the structure into a 
finite number of small elements. Within each element, the dis- 
placements u and v at each point (x,y) in the x- and y-direction 


are assumed to be series functions of x and y of the form: 


i ae 
- 2 1 J : : 
u(X,y) é ; as Key 
v(x,y) = 2 2b.. xt y! <n _ 
(x,y) he y Lo. 
Bas and Dat are functions of the coordinates (py Yy) of the nodes 
My limiting the elements. The functions u(x,y) and v(x,y) are such 


that displacement compatibility is satisfied along the boundaries 
of the elements to ensure the condition of convergence to the 
true solution. By the use of the variational method, a system 


of linear equations between forces and displacements can be derived 


l 





® and can be solved for forces if displacements are given or vice 


Velsda. 


The computer programs set up by the Civil Engineering Depart- 


ment for plates are of four types: 


(1) "CSTG": constant strain triangle, global formulation in which 
the elements are triangles, the strain within each element is 
assumed to be constant and one global referential system is 
used for the whole plate. 

(2) "CSTL": constant Strain triangle, Local Lormulat1onaeeeie 
only difference with the "CSTG" is that each element has its 
Own individual referential system. 

(3) “LST": linear strain triangle. The elements are triangles, 


the strain within each element is assumed to vary linearly. 





(4) "PSR": plane stress or plane strain rectangle tice te wes 


are rectangles. 


The finite element methods provide a powerful means of anal- 
ysis of structutte. Any complicated structure can be investigated 
without much difficulty. They have, however, their own weaknesses. 
The input data need to be carefully set up. The results are ob- 
tained in the form of numerical answers to the problem and the 
influence of different Eactors are not evidenced liners toner 


where it is possible, analytical approach will be used. 


II Choice of Appropriate Numerical Method 


There are two alternatives in the use of the previously- 

















mentioned computer programs, namely the constant strain ele- 


ment: "CSTG", "CSTL", “PSR” or the linear strain. eleneninus 5 oe 


Therefore the computation of the degree of constraint in 
the case of a straight slit in a rectangular plate has been 
made first, using both the constant strain element and the linear 
strain element programs as means of verifying tne adacmacs Ou 
the methods themselves and a comparison between the two program 
byes Since analytical results are available (equation[3], for 
an infinite plate). Also in the case of constant strain, two 
types of gridwork are used: one relatively coarse, made of a 
combination of triangles and rectangles, CSG “and sh eey oor 
with 105 nodes. The other 1S made with a finer gridwork, of 
triangles, "CSTG" type, with 170 nodes. The results are pre- 
sented in Tables 1 to 6, the loading being symmetrical with 
respect to the center of the slit. The input and output data 


are presented in Appendices la and lb for constant strain and 


in Appendix 2 for linear strain. The Symbols: used are: 


R= 


al ho 


£: loading Length 
L: length of the slit 


v: local transverse displacement (in the direction 
perpendicular to the slit) 


average transverse displacement over the loading 
length & 


K : degree of constraint computed 


K = —<e : 
Iv}, 








. 
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K : non-dimensionalized degree of constraint computed: 


K E: Young Modulus of 
[ee Elasticity 
E/L L: Length of the slit 


K_: non-dimensionalized degree of constraint obtained 
analytically in the case of an infinite plate: 


The use of non-dimensionalized K (K. or K.) makes the results 
somehow comparable in the form of curves relating K or Ke to R and 
are presented in Figure 5. The difference between K, and K could 
not be considered as a measure of the accuracy of the numerical 
method because Ky is the non-dimensionalized degree of constraint 
in an infinite plate and for which the only significant dimension 
is the slit length Lb of the welds ances ae non-dimensionalized 
degree of constraint ina finite plate and for which the signifi- 
cant dimensions include the slit length L, he plate length L, 
and the plate width W,. But we can expect K to be smaller than 


K, since the restraint due to a finite plate sheula@ bc ole wera 
; 


that due to an infinite plate. 
The average transverse «displacement [vl], is defined-as: 


L/ 2 
vl, =7 f vax 
yO 


for a symmetric weld, where 2 is the weld length. 


In Onder stoudete anuinieomi the local waluesmoty [ome 


Ll 


ious nodes where loads are applied,are plotted as 4 £unctrvonges 








the abscisse x o£ the nodes. Then a Curve 16) stawmte cea 
through these points and the area limited by the axis Oy 

(which is also the axis of symmetry), the ordinate x = 8/2, 
the curve (C) and the axis 0x is measured graphically by a 


planimeter. 


Area= vax 


zee 


A,B,C,D are loaded nodes 
Me 
D 2 
By dividing the value of this area to 2/2, the average 
half transverse displacement v/2 is obtained, since the weld 
is also symmetrical with respect to 0x and the transverse 
displacement is the distance which both sides of the weld 


come in to each other. 


From the results in Figure 5 these conclusions could 


be formulated: 


(1) The linear type gives an unusual swing: the curve of K 
is above that of K. for small R and becomes smaller than 
that of K computed by the "CSTG" method for large va- 
IMeSyvOren. 

(2) The constant strain types "CSTG", "PSR" in the two 
gridworks give more consistent results. The K curve 
obtained is almost parallel to the K, curve. Also, 


there is not much difference between the two gridworks, 
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One COarsSe ona thermo tnenaem= anes 


Therefore, the constant strain element type with a 
relatively coarse gridwork can be used in the computation 
of the degree of constraint of various joint configurations 


in« plate Stpuceume. 


III Units Used 


1) Length: 


(a) inch 
(6b) millimeter wherever comparison with Japanese data 


1s desirable. 
2) Otce. 


(a) Pound 
(b) Kilogram wherever comparison with Japanese data 


1S Cle eSics ale - 
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Dimension Chatacterisel cs. 
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ie . ' = b , 3 7 

Length of Plate,L,=6" Length of Sit, )i=2: 

Width of Plate,W,=4" Maximum Width of S1lit,W=0.2" 
Plate Thickness ,T=0.2" web SSE LS 10 


Small Axis b 


Mechanical Characteristics: 


Force Applied. — (2eec) 10. ello, ann 


O Oe socul 


Yo 


Young's Modulus, y=" 50)bal0 ect 
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Values of the Degree of Constraint K and the Non-Dimensionalized 


Degree Of GonctEaen aan 


—— 
(E/L) 


TABLE 2 


(Constant Strain Element) 


v/2 


TMOG sin) 


1.797 
95128 
4.192 
5.039 
5.781 
6.315 
6.730 
6.904 


6.866 


K 


(10= Sosay isa) 
27.824 
16.000 


LSS 2 


Seize 
8.649 
Toe lke 
T2429 
Tea ae 


RPS? 


AI 


Lo BSS 
1.067 
0.795 
amGe? 
VST 
0.528 
0.495 
0.483 


0.485 


for Different Values of R 


A | 


2 O99 
Ve ZOS 
Wid 2 
Oo23 
Os 7) Sil. 
oy 2 
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Oo iy 
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TABLE 4 


Values of the Degree of Constraint K and the Non-Dimensionalized 





Degree of Consteraln oe e— TEEY for Different Values of R= 2 
(Constant Strain Elements - Finer Gridwork) 
9 v/2 K _ _ 
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Oa 125 28.986 ES) 2 | a Oe, 
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RESULTS 


L Presentation of the Results 


A. Straight Slit Type - The results are assembled in tables 
7 and 8 and in Figure 6; the input and output data are in 


Appendix 3. 


B. Straight Slit with Circular Holes at the Ends —- The re- 
Sults are in Tables 9 and 10 for the "CSTG" type and in 
Tables 11 and 12 for the "LST" type and are assembled in 
Figure 7 for both types. The input and output data are in 
Appendices 4 for "CSTG" and 5 for "LST." Here again the 
two methods, "CSTG" and "LST" are used in order to compare 
them again because the weld configuration is different from 


the previous cases. 


C. H-Slit Type, Steel Plate - The results are in Tables 13 
and 14 and ,n Figure 8.9)7he input and eutpuT data a ee 


Appendix 6. 


D. Lehigh Test Specimen - The results are in Tables 15 and 


16. The input and output data are in Appendices 7a and 7b. 


E. H-Slit Type, Aluminum Plate - The ~esults are inelab lee 


The input and output data are in Appendix 8. 
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TABLE 8 


Values of the Degree of Constraint K and the 


Non-Dimensionalized Degree of Constraint K = TELY 


for Different Values of R = 2 
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TABLE 10 


Values of the Degree of Constraint K, the 


Non-Dimensionalized Degree of Constraint K = a. 
and the Specific Degree of Constraint K' = - 
for Different Values of R = - 
(Constant Strain Elements) 
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Node Displacements in the Y-Direction (10. 
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TABLE 14 


Values of the Degree of Constraint, the 


Non-Dimensionalized Degree of Constraint K = SIT 
and the Specific Degree of Constraint K' = K/E 
for Different Values of R = &/LL 
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Length of Plate, L,=300 mm Slit Overall Length,L=150mm 
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Plate Thickness, T=10 mm End Circle Diameter ,D=12.5mm 
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TABLE 17 


Result of the H-Slit, Aluminum Plate 
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II. Discussion Of Results 


ae Interpretation of the Results 


The valid results for different joint configurations 


are assembled in Figure 9, 


As expected, the K value for a finite plate is smaller 
than that for an infinite plate. The effect of finiteness 
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Both plates have the same aspect ratio L,/W, = 3/2. It can 
therefore be concluded that for straight slit weld config- 
uration, (there is not much difference between an ellipse 
with an axis ratio of a/b = 10 and a straight Slit ene 
value of K depends on the ratio between plate length L, and 
Slit length L. K is larger for larger ratio L,/L and for 
straight slits having the same length, the degree of con- 
straint K is larger for the one with larger 7 len eae 
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In the case of a straight slit with circles at the ends, 
the diameter of the circle is larger than the slit width, ene 
value of K is smaller than that for a straight slit without 
large-circled ends. And the lowest value of K is obtained 
for an H-slit, where the width of the H-branch is larger than 


that of the weld slieybW — 107 3>- 2-2 


All these results seem to conform to common sense and 
could be considered as an indication of the correctness of 
the procedure. Theoretically, it cannot be otherwise since 
it can be proved that, in the finite element method, if the 
compatibility conditions are satisfied, the results will 
converge to the true solution as the gridwerk becomeo ei. = 
and finer. In the cases of constant strain element, the 


compatibility conditions are satisfied. 


A full verification of the method is made in the case of | 
the Lehigh specimen (tables 15 and 16). In the case "without 
Mehifete ils eeewlics! eugee 


= 48.17 Kg/mm-mm , = 44 Kg /mm-mm 


AK _ 4.17 


K 48.17 


K K 
measured 


7.17 te) Inv ehe Gases ise 


calculatéd 


and the difference is 
sawcut", they are: 


= 30.69 Kg /mm*—mm , = 27 Kg /mm*—mm 


MORE Ves 
K 30.69 


RAS Irene 
2 O22 


Kalculated 


II 


ana eherarrtrerence mms 


These values of K have been computed by taking the 


3 


Young's Modulus E = 21.1 x 10 Kg/mm", which is a little high. 


| 
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If we take E = 19.9 x 10° Kg/mm“, then we will obtain 


ee 2 leg, 
K' = a1. 1 ES) 2S, 





45.3 Kg/mm--mm 


for the plate without cut and 


xt = 30:69 x 19.5 





28.9 Kg /mm—mm for the plate with sawcuf. 


With these values, the differences between the calculated and 


measured values of K become: 


ORL _ 4523-44 _ 1.3 _ 5 ge 

Ky a5.3 45,3 > 7°97 
and 

e2 _ 28.9-27 | 1.9 pone 

K, 28.94 ~— 28.9 sal 


With these differences, the values of K calculated seem 
to be acceptable if we think of the complicated boundary con- 
‘dition in the case of the plate with sawcut as described on 


page 60. 


B. Significance of the Results 


Coming back to Figure 9, it can be seen that the value 
K increases sharply for a value of x = £/L smaller than 
about 0.30. The degree of constraint K is related to the 
transverse stress Sue by . Kiv],- From the empirical 
Watanabe-Satoh relation, the transverse shrinkage is related 


Oe 2G Sh 


7 
SSG) ie 
C tr i+0.086 Ko°8! 
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Along the weld line Sa = v and for a long straight weld, 





val ~ vy since v is almost uniform, except near the ends, 
then, by substituting the expression of Sy in the relation 


giving Oye we have: 
oO 


— ——~ K 
5 = See See 
Yo tr i+0.086 KO??? 


where o,, 1s the average reaction stress along the weld due 
‘ © 


to constraint. 


The value of O., /See as a function of K is plotted below. 
oO 


——s 
—- 
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oO 
Yo 
dered as an indication of the residual stress level near the 


iS an increasing function Gf K, and couleebeoecencn. 


weld. 


Therefore, it Gan be said that 1£ aust rerol amos eenee 
large plate is to be welded, and if the weld is long enough 
so that more than one block is necessary, the residual stress 
would be lower if each weld is made symmetrically and if the 
ratio of the weld length 2 to the slit length L is larger 


Ginteamy ADOU te asl oe 


Also, the Watanabe-Satoh relation can present a problem 
of choice: if the transverse shrinkage is to be limited to 
small values, then high eee of constraint and consequently 
high level of residual stress will result. On the contrary, 
if the residual stress is to be small, then low degree of 
constraint 1s necessary (if the degree at constraint can be 
controlled) and high values of transverse shrinkage will 


result. 
{ 4 


On the other hand, welded structures are subjected to 
cracking. Mechanical factor is one of the main factors nae 
promote cracking. The reason might be attributed to high 
residual stress associated with high degree of restraint. 
Therefore, it may be important to determine a critical value 
of the degree of restraint K. necessary to produce cracking. 
Experimentation iS necessary and a series of tests, known as 


the Lehigh test specimen (reference 4, page 32) is proposed 
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9 to determine KL: Basically, it consists of several plates 


with a straight weld made along a part of its longitudinal 


median line. If the plate is not cut, a high degree of con- 
| | j . <= 
: ( a) | Ces i | 
| | f | ' 
Soke _ |_| | te 


straint can be achieved and cracking may occur. To decrease 
the degree of constraint of the weld, several sawcuts are 
made along the perimeter of the plate, and the degree of con- 
Straint decreases as the length of the sawcut increases. 

Then there should be ee ee of the sawcut length for 

, which cracking begins to disappear, the corresponding value 
of the degree of constraint 1s the critical value K, which 


can be considered as a material constraint for a given weld 
iy Cr 

Pius Wikeho the value or K, determined, we can predict 
whether a weld will crack or not due to its mechanical re- 
straint, Since the value of K can be computed for a given 
weld configuration. One field of application could be in 
the patch weld, which consists of welding a circular disc 


= ae aa 


3 r | 
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in a flat plate. If the dimension of the hole is such that 
the disc dimension is small compared to the plate dimension 
and the weld can be made in one pass, then we have (reference: 
Mathematical Theory of Elasticity, by Sokolnikoff, Chapter 5): 


20,°R 
Sp = 3 (for an infinite plate) 





Op is the uniform radial stress at the weld, 


Sp is the radial shrinkage along the weld and is 
equal to the sum of the radial displacements of the disc and 


of the plate at the radial distance, r = R. 


Then: 


an 

lI 

lI 
No] 
73 | bd 


In this formula, we find that the degree of constraint 
in the case of a complete patch weld is inversely proportional 
to R. Therefore, there should be a critical value Ry Of kK 


below which K > K, and the weld will crack. 
However, this value of K: 


—— 


Se 
ee) 


is an upper limit, for it is exact only for an infinite plate. 
For a finite plate, the value of K will be small. The effect 
of relative dimension could be seen more clearly by examining 


the expression for the radial stress outside the disc, which is: 
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At infinite, oe Oe sLE R/ raed 710,. then J /Op =<) / A O0vee te 
a distance equal to 10 times the radius of the disc, the 
radial stress is one hundredth of the stress at the weld. 


This value might be still important in some cases. 


If the disc is large so that the weld can only be made 
sector by sector, and if the weld is made symmetrically with 


respect to the disc center, then with an infinite plate, we 





é 


[5 nt A OA 


| 


- 
u 
. 


have, for-a weld sector of angular value 20,as in the figure: 


O 

2 R LE il 

K (20) = — = ik Aus oe 
[up] 


20 
where £(20) 1S a convergent series (it converges as Wines with 


© (1h) oo — ie (0) On ee is the mean value of the radius dis- 


IAC! 


placement taken along the weld angular distance of 2a. 


UR 


The form of the expression suggests that for K(2a) to be 


smaller than a critical value Kt 
< 
K (2a) K, 
then for a given radius R, there is a value 20. EO ewhie in 
K (2a) = Ky 


and for a < a, KK. Thus, if a large circular disc is to 
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be welded to a plate and the weld has to be made sector by 
sector, symmetrical around the origin, then there is a cri- 
tical value 20, of the angular weld sector below which the 


weld might crack and above which the weld would not crack. 


III. The Problem Of Experimentation 


In the past, experiments have been conducted by actually 
welding the specimens and measuring the transverse shrinkage 
and the value of the degree of constraint. Little use has 
been made of the elastic assumption and of the resulting 


linear relations between stress and displacement. 


One inconvenience of the measurement of the strain due 
to welding is the high temperature at the weld, and strain 
gauges had to be put some distance away from the weld. On 
the other hand, because the weld line is not well defined 
after the weld is made, the displacements had to be measured 
also some distance away from the weld. Furthermore, in some 
cases, weld specimens are difficult to make, as the case of 


a small slit ina big plate,and consequently expensive. 


Taking into account the elastic assumption, it can be 
noted that for a given average transverse stress, Or at the 
weld line, the local displacement and hence the average dis- 
placement is inversely proportional to the Young's Modulus E. 
Therefore, the degree of constraint, defined as: 

[vl], 








Ue 


is proportional to E; ‘and the ratio of K and eb, whieh ia. 





been defined as specific degree of constraint K': 


K 
oe 
" E 

is independent of the material, provided that it is linear 


elastic. 


Consequently, there could be a simple method of experi- 
mentation in the measurement of the degree of constraint K by 
having specimens made of material that has a low modulus of 
elasticity E and by not actually performing the weld but by 


using some mechanical means to simulate the thermal stress. 


One such material could be plastic such as cast phenolic 
% resins for mechanical and chemical purpose whose stress- 
strain curves are reproduced in Figure 10 (reference: "Tech- 
nical Data on Plastics} by Macias £ Selena Chemists' Associa- 
tion, Inc., page 49) and which have a modulus of elasticity 


in compression of 3 x 10> psi to 5 x 10> psi (ibidem, page 52). 
; 


To impose stress in the model, wedges can be driven into 


the slit. In this case, two similar prismatic wedges with 
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elongated rectangle triangle cross section made by cutting 

a prismatic wedge with rectangular cross section along the 
Giagonal plane should be used so that the faces in contact 
with the slit edges remain parallel. The sample will be in 
compression to minimize the friction between the active faces 
of the wedges and the slit edges. Some lubricant, such as 


Teflon, graphite or oil could be used. 


a ee ~ 
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ft 
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at cea dees 
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The width of the wedge in contact with the plate repre- 


sents the weld length &. Several sets of wedges with dif- 

ferent widths are necessary to establish a curve giving the 
specific degree of constraint K' as a function of the ratio 
R between weld length 2 and slit length L: R = &£/L for each 


weld configuration. 


{ s 
To measure the stress, sensible strain gauges can be 


used along the edges of the slit on both sides. To measure 
the transverse displacement along the slit, visible marks 
are made along the edges of the slit and their displacements 
can be measured by a microscope, or a comparator. After 
measurement of local stresses and displacements along the 
load length &, their averages will be computed and the spe- 


cific degree of constraint K' obtained by :; 
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For the verification of the Watanabe-Satoh curve, there 
is,however, unfortunately no simple way and full-size experi- 
ments have to be performed. Given a Specimen, we can first 
compute its degree of constraint and deduce its transverse 
shrinkage uSing the Watanabe-Satoh curve. Then actual wel- 
ding is performed and experimental data obtained are com- 


pared with computing results. 


If the Watanabe-Satoh relation is correct, then the 
knowledge of the value of K is most important in the trans- 
verse shrinkage control Since these laboratory results seem 
” to be applicable to actual structures as has been mentioned 
in reference 3. In these Situations, if the degree of con- 
Straint can be computed, then the transverse shrinkage can 
be determined and hence the level of residual stress is also 
known which will enable one to predict whether weld cracking 
due to mechanical constraint could happen and subsequently 


to take necessary corrective steps if possible. 














CONCLUSION 


The finite element method, with constant strain ele- 
ments, seems to give consistent results. A qualitative 
verification of its validity remains somehow in the fact 
that the results conform to common sense. Two quantitative 
verifications have been made and the results can be consi- 


dered to be acceptable, in the case of the Lehigh specimens. 


Taking into account the elastic assumption, and the 
resulting proportionality between the degree of constraint 
K and the modulus of elasticity E, a method of experimen- 
tation that enables the determination of k experimentally 


has been proposed. 


The determination of the transverse shrinkage due to 
butt weld of a plate structure can be made, if the degree 
of constraint is known, by using the Watanabe-Satoh corre- 


lation between relative transverse shrinkage S/S and de- 


feat 
qi CC Of Be@oms iad aie 

The Watanabe-Satoh correlation has been obtained empiri- 
cally through experiments performed in three different spe- 
cimens: the straight slit, the H-slit and the ClrCUulargisimae 
some further experiments on various slit types would be 


desirable. 


If the numerical method proves to be good and the 


Watanabe~Satoh correlation to be valid, then shrinkage dis- 











te, 


tortion control 1s possible. In this case, theres ateuc ie 
alternatives, deduced from the Watanabe-Satoh curve: either 
to accept a high degree of constraint and hence a high level 
of residual stress and a low value of transverse shrinkage 
or a low value of the degree of constraint and hence a low 
level of residual stress and a high value of transverse 


shrinkage. 


On the other hand, welded structures are subjected to 
cracking. External constraint could be one major mechanical 
factor contributing to weld cracking. For each material and 
weld type, there could be a critical degree of constraint, Kae 
above which weld cracking may occur and below which it may 
OKO, - mhe Lenton test specimens are used for this purpose. 
Once the critical value K. is determined, it will be possible 
to Know in advance whether a weld is crack-susceptible and 
therefore to decide on its feasibility or corrections if 
necessary and possible. One example of application can be 


rd 


found in patch welding. 


In both cases, distortion control or weld cracking sus- 
ceptibility verification, it is essential to compute the 


value of the degree of constraint K of the weld. 
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